siml.c 

#include <stdio.h> 
int 

main() 
{ 

int k; 

double pi , p2 , p3 , p4 , p5 , p6; 
double al , a2; 
double ml , m2; 

double const CONVERGE_TOL =0.1; 
double delta , old_val; 
int i_count=0; 

pi = 0; p2 = 0; p3 =0; p4 = 0; p5 =0; p6 = 0; 
delta = lel6; 

while( delta > converge_tol ) 
{ 

old_val = pl-p2; 

/* apply rule "I" */ 
p2 -= 1000.0; 
pi += 1000.0; 

/* apply rule "x" to intersection cells */ 

al = (pi + p3 + p5)/3; 

pi = p3 = p5 = al; 

a2 = (p2 + p4 + p6)/3; 

p2 = p4 = p6 = a2; 

/* apply rule "r" to resistive elements; */ 

/* specifically Rl = 200 Ohms, R2 = 400 Ohms */ 

ml = (p3 - p4)/200.0; 

p3 -= ml; 

p4 += ml; 

m2 = (p5 - p6)/400.0; 
p5 -= m2; 
p6 += m2; 

fprintfC stdout , "pi %10.2f\tp2 %10.2f\tpl-p2 %10.2f\n" , pi , p2 , pl-p2 ); 
fprintf( stdout , "p3 %10.2f\tp4 %10.2f\tp3-p4 %10.2f\n" , p3 , p4 , p3-p4 ); 
fprintf( stdout , "p5 %10.2f\tp6 %10.2f\tp5-p6 %10.2f\n" , p5 , p6 , p5-p6 ) ; 

/* check for convergence */ 
delta = (pl-p2) - ol d_val ; 
i_count++; 

fprintfC stdout , n iteration[%d] delta=%g \n\n" , i_count , delta ); 



} 



} 



/* 

pi 333.33 p2 -333.33 pl-p2 666.67 

p3 330.00 p4 -330.00 p3-p4 660.00 

p5 331.67 p6 -331.67 p5-p6 663.33 
iteration[l] delta=666.667 

pi 665.00 p2 -665.00 pl-p2 1330.00 

p3 658.35 p4 -658.35 p3-p4 1316.70 

p5 661.67 p6 -661.67 p5-p6 1323.35 
iteration[2] delta=663.333 

pi 995.01 p2 -995.01 pl-p2 1990.02 

p3 985.06 p4 -985.06 p3-p4 1970.12 
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p5 990.03 p6 -990.03 p5-p6 
iteration [3] del ta=660 . 017 
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pi 66656.69 p2 -66656.69 pl-p2 

p3 65990.12 p4 -65990.12 p3-p4 

p5 66323.40 p6 -66323.40 p5-p6 

iteration [1757] del ta=0. 100286 



pi 66656.74 
p3 65990.17 
p5 66323.45 
iteration [1758] 
*/ 



p2 -66656.74 pl-p2 

p4 -65990.17 p3-p4 

p6 -66323.45 p5-p6 
del ta=0. 0997842 



133313.38 
131980.24 
132646.81 



133313.48 
131980.34 
132646.91 
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